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Superpositioning  of  relaxation  data  as  a  function  of  the  product  variable  TV7 ,  where  T  is  temper¬ 
ature,  V  the  specific  volume,  and  y  a  material  constant,  is  an  experimental  fact  demonstrated  for 
approximately  100  liquids  and  polymers.  Such  scaling  behavior  would  result  from  the  intermolecular 
potential  having  the  form  of  an  inverse  power  law  (IPL),  suggesting  that  an  IPL  is  a  good  approxi¬ 
mation  for  certain  relaxation  properties  over  the  relevant  range  of  intermolecular  distances.  However, 
the  derivation  of  the  scaling  property  of  an  IPL  liquid  is  based  on  reduced  quantities,  for  example,  the 
reduced  relaxation  time  equal  to  7'//2V  ' 1/3  times  the  actual  relaxation  time.  The  difference  between 
scaling  using  reduced  rather  than  unreduced  units  is  negligible  in  the  supercooled  regime;  however, 
at  higher  temperature  the  difference  can  be  substantial,  accounting  for  the  purported  breakdown  of 
the  scaling  and  giving  rise  to  different  values  of  the  scaling  exponent.  Only  the  y  obtained  using 
reduced  quantities  can  be  sensibly  related  to  the  intermolecular  potential.  ©  2011  American  Institute 
of  Physics,  [doi:  10. 1063/1.3532545] 


Extensive  measurements  over  the  last  few  years  of  re¬ 
laxation  properties  of  organic  materials  at  elevated  pressures 
have  led  to  many  insights  into  the  dynamics  of  glass-forming 
liquids  and  polymers. 1,2  An  interesting  aspect  of  these  studies 
is  the  observation  that  over  a  wide  range  of  thermodynamic 
conditions  relaxation  times  conform  to  the  scaling  law 

r  =  f(TVn,  (1) 


where  /  is  a  function,  V  the  specific  volume,  and  y  a  mate¬ 
rial  constant. 3-6  Equation  (1)  has  been  extended  to  other  dy¬ 
namical  quantities,  such  as  the  viscosity  t]1,&  and  diffusion 
coefficient  D9  and  also  applied  to  molecular  dynamics  sim¬ 
ulations  (mds)10  15  of  particles  interacting  via  Lennard-Jones 
(LJ)-type  potentials 


U(r)  =  e 


(7)" -  (7/ 


(2) 


where  0  is  the  particle  diameter,  e/4  the  well  depth,  and  n 
a  material  constant.  Of  particular  interest,  mds  results  have 
demonstrated  a  connection  between  the  y  in  Eq.  (1)  and  the 
steepness  of  the  intermolecular  repulsive  potential  [related  to 
the  n  in  the  repulsive  term  in  Eq.  (2)]. 12,15  This  is  intrigu¬ 
ing  because  the  idea  underlying  density  scaling  is  that  the  dy¬ 
namic  properties  of  viscous  liquids  are  governed  primarily  by 
the  repulsive  component  of  the  intermolecular  potential,  with 
changes  in  the  long-range  attractive  part  exerting  a  negligi¬ 
ble  effect.  To  the  extent  this  is  the  case,  for  certain  properties 
Eq.  (2)  can  be  approximated  by  an  inverse  power  law  (IPL) 

U(r )  =  £  (-)  .  (3) 

Since  Eqs.  (2)  and  (3)  are  functions  of  distance  but  not  ori¬ 
entation,  such  potentials  can  only  be  approximate  for  non- 
spherical  molecules.  However,  the  latter  can  be  described 
using  multiple,  discrete  LJ-site  interactions, 16-1 8  and  such 
models  exhibit  correlation  of  fluctuations  in  potential  energy 


and  the  virial.14  This  is  consistent  with  an  effective  IPL,  since 
the  potential  energy  and  the  virial  are  perfectly  correlated  only 
for  an  IPL.  This  applicability  of  the  IPL  approximation  to  the 
supercooled  dynamics,  however,  does  not  imply  any  general 
validity.  For  example,  the  melting  behavior  of  nonspherical 
molecules  exhibits  strong  deviations  from  density  scaling. 19 

The  excess  transport  and  thermodynamic  properties  (“ex¬ 
cess”  refers  to  the  nonideal  part,  dependent  on  particle  po¬ 
sitions  but  not  momenta)  of  an  IPL  liquid  are  a  function  of 
TV"/ 3. 20,21  -phis  follows  from  the  fact  that  an  IPL  and  its 
derivatives  depend  only  on  the  combined  variable  £er~",  so 
that  properties  depend  only  on  a  single  temperature-density 
variable.  The  implication  is  that  conformance  of  real  materi¬ 
als  to  density  scaling  [Eq.  (1)]  results  from  their  intermolec¬ 
ular  potential  being  approximately  an  IPL,  at  least  for  con¬ 
sideration  of  certain  “local”  dynamic  properties.  Dyre  and 
coworkers14, 15  coined  the  term  “isomorphic  states”  to  refer  to 
state  points  in  which  the  partition  function  constructed  from  a 
given  potential  depends  to  a  good  approximation  (or  exactly 
for  an  IPL)  on  the  single  variable  TV' . 

Although  for  a  strict  IPL  y  =  h/3,  results  for  LJ  liquids 
indicate  that  the  scaling  exponent  is  larger  than  the  repulsive 
exponent  (y  >  n /3),12  the  steeper  repulsive  slope  (h  >  n)  due 
to  the  presence  of  the  attractive  term  in  Eq.  (2). 22  Neverthe¬ 
less,  such  results  suggest  that  density  scaling  of  experimental 
data  has  the  potential  to  yield  direct  information  about  the 
forces  between  molecules  in  real  materials.  However,  there  is 
a  disconnect  between  scaling  of  dynamic  data  from  mds  and 
from  experiments.  For  the  latter,  scaling  is  applied  to  r,  r], 
and  D  as  measured,  with  the  scaling  property  being  an  ob¬ 
servable  fact.  However,  the  scaling  property  of  effective  IPL 
potentials  applies  to  reduced  quantities,  using  units  of  time 
to  =  v1/3  ( kT /m)~l/2.  length  l0  =  v1/3,  and  energy  E0  =  kT, 
where  m  and  v  are  molecular  mass  and  volume.20  21  Dimen¬ 
sionless  reduced  quantities  are  thus  defined  as 
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x*  =  v  l^(kT /m)l^2r, 
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rf  —  v2^(mkT)  l^2ri, 

D*  =  v~lp(kT/m)-1,2D.  (4) 

Accordingly,  for  an  IPL  and  mds  data  the  form  of  Eq.  (1)  is 

r*  =  f(TVy\  (5) 

with  similar  expressions  for  rj*  and  D*.  Since  over  the 
range  of  experimental  measurements  on  supercooled  liquids, 
T  varies  by  as  much  as  100%  and  v  by  ca.  10%,  it  is  not 
a  priori  obvious  that  differences  between  reduced  and  nonre- 
duced  units  are  negligible.  Conformance  of  experimental  data 
to  Eq.  (5)  and  any  difference  between  y  *  and  y  are  important 
issues  if  the  scaling  exponent  is  to  be  interpreted  as  the  slope 
of  an  effective  intermolecular  potential. 

Density  scaling  is  carried  out  by  plotting  a  dynamic  prop¬ 
erty  variable  versus  T\r/ ,  with  y  adjusted  to  collapse  the  data 
onto  a  master  curve.  If  the  variation  in  r,  rj,  or  D  is  much 
greater  than  that  of  volume  and  temperature,  then  the  shape 
of  this  master  curve  will  be  essentially  identical  whether  ac¬ 
tual  or  reduced  quantities  are  employed.  Since  the  hallmark 
of  the  supercooled  state  is  the  enormous  change  in  molecular 
mobilities  for  small  changes  in  thermodynamic  conditions,2 
scaling  of  r  or  r*  is  practically  equivalent.  This  is  illustrated 
for  squalane  in  Fig.  1,  which  is  representative  of  the  literature 
results  for  dielectric  relaxation  times  of  supercooled  liquids 
and  polymers. 

Many  liquids  cannot  be  supercooled  and  relaxation  data 
(typically  viscosity  or  self-diffusion  coefficient  measured  us¬ 
ing  NMR)  are  available  over  only  a  narrow  range  above  the 
melting  point.  In  such  circumstances  the  use  of  reduced  units 
makes  a  substantial  difference.  Fig.  2  shows  high-pressure 
viscosity  data  for  a  typical  case,  dodecane.  Pensado  et  alf 
reported  a  scaling  exponent  of  y  =  6.5  for  the  viscosity, 
with  good  superpositioning  for  smaller  values  of  TV y,  but 
the  scaling  deteriorated  for  lower  rj.  Adjusting  the  value  of 
y  improves  scaling  for  lower  rj,  but  then  the  higher  viscosity 
data  diverge.  Alternatively,  the  scaled  viscosities  are  plotted  in 
Fig.  2,  and  the  superposition  is  very  good  over  the  entire  range 
with  y*  =  5.2.  Note  that  if  the  dynamic  data  encompassed  a 
narrower  range,  as  is  often  the  case,  both  reduced  and  nonre- 


FIG.  2.  Density  scaling  of  viscosity  (open  symbols)  and  reduced  viscosity 
(filled  symbols)  of  dodecane.  The  reduced  viscosity  shows  better  scaling  and 
a  lower  scaling  exponent. 


duced  viscosities  would  ostensibly  scale  equally  well,  albeit 
with  different  values  of  the  scaling  exponent.  For  example,  for 
toluene  y  —  7. 8, 7  whereas  we  find  y*  =  5.4;  for  n-hexane  y 
=  13, 8  while  we  find  y*  =  1.  Also,  the  scaling  exponent  y 
is  independent  of  the  experimental  variable  (relaxation  time, 
transport  coefficient,  etc.)  only  when  reduced  units  are  used. 
This  is  illustrated  in  Fig.  3  for  the  viscosity  and  diffusion  con¬ 
stant  of  cyclohexane  (using  data  from  Ref.  23). 

We  can  also  apply  the  density  scaling  relation  to  very  low 
viscosity  liquids,  such  as  krypton  (shown  in  Fig.  4),  argon, 
nitrogen,  oxygen,  and  methane  (using  viscosity  data  from 
Ref.  24).  For  these  substances,  nonreduced  viscosities  do  not 
superpose  versus  TVy  for  any  value  of  y.  However,  the  re¬ 
duced  viscosities  conform  to  the  scaling  relation,  with  physi¬ 
cally  reasonable  values  of  y  *  =  4.0  for  krypton,  argon,  nitro¬ 
gen,  and  methane,  and  5.0  for  oxygen. 

The  above  results  affirm  the  connection  between  y  and 
the  steepness  of  the  intermolecular  potential  for  r  relevant 
to  the  local  dynamics.  Density  scaling  is  not  merely  an  em¬ 
pirical  relation  useful  for  organizing  experimental  data,  but 
has  a  physical  basis  in  the  IPF  approximation  for  the  in¬ 
termolecular  potential.  If  the  data  encompass  a  range  suf¬ 
ficient  to  yield  different  values  for  y  and  y  *,  the  latter, 
based  on  reduced  dynamic  quantities,  is  the  physically  rele¬ 
vant  exponent.  Examples  for  which  y  *  differs  from  published 
y  include  n-alkanes,8  octane,  and  toluene,7  1-alkylamines,25 


FIG.  1.  Density  scaling  of  viscosity  (open  symbols)  and  reduced  viscosity 
(filled  symbols)  of  squalane.  The  data  extend  over  many  decades  of  viscosity, 
therefore  both  quantities  scale  with  identical  scaling  exponents.  In  this  and  all 
other  figures  TV7  is  in  units  of  K  (ml /g)Y  and  viscosities  are  mPa  s;  reduced 
units  are  dimensionless. 


FIG.  3.  Density  scaling  of  the  viscosity  and  diffusion  coefficient  of  cyclohex¬ 
ane  using  unreduced  (open  symbols,  D  in  cm2/s)  and  reduced  (filled  symbols) 
units.  For  the  latter  the  exponent  is  independent  of  the  dynamic  variable. 
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FIG.  4.  Density  scaling  of  viscosity  (open  symbols)  and  reduced  viscosity 
(filled  symbols)  of  krypton.  The  reduced  viscosities  scale  well  with  an  expo¬ 
nent  y*  =  4.0;  however,  the  scaling  is  poor  for  any  value  of  y . 

and  2-alkylamines,26  carbon  dioxide,  and  toluene.9  Moreover, 
some  examples  of  an  ostensible  breakdown  of  the  scaling  are 
a  consequence  of  using  nonreduced  variables. 
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